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1. Introduction 

Experiments show that many excited-state hadrons exist, and there are significant experimental 
efforts to map out the QCD resonance spectrum, such as Hall B and the proposed Hall D at Jefferson 
Lab, ELS A associated with the University of Bonn, COMPASS at CERN, PANDA at GSL and 
BESIII in Beijing. Hence, there is a great need for ab initio determinations of such states in lattice 
QCD. 

Higher-lying excited hadrons are a new frontier in lattice QCD, and explorations of new fron- 
tiers are usually fraught with dangers. Excited states are more difficult to extract in Monte Carlo 
calculations; correlation matrices are needed and operators with very good overlaps onto the states 
of interest are crucial. To study a particular state of interest, all states lying below that state must 
first be extracted, and as the pion gets lighter in lattice QCD simulations, more and more multi- 
hadron states will lie below the excited resonances. To reliably extract these multi-hadron states, 
multi-hadron operators made from constituent hadron operators with well-defined relative mo- 
menta will most likely be needed, and the computation of temporal correlation functions involving 
such operators will require the use of all-to-all quark propagators. The evaluation of disconnected 
diagrams will ultimately be required. Perhaps most worrisome, most excited hadrons are unsta- 
ble (resonances), so the results obtained for finite-box stationary-state energies must be interpreted 
carefully. 

This talk will describe the key issues and challenges in exploring excited hadrons in lattice 
QCD, emphasizing the importance of multi-hadron operators and the need for all-to-all quark 
propagators. Dealing with resonances in a box is discussed, and the technology associated with 
extracting excited stationary-state energies, including operator design and field smearing, is de- 
tailed. Efforts in variance reduction of stochastically-estimated all-to-all quark propagators using 
source dilutions are outlined. Results on excited hadrons during the last year are summarized. 

2. Resonances in a box 

A simple example serves to illustrate the issues that must be confronted when studying reso- 
nances in a box. Consider the (dimensionless) Hamiltonian for a single nonrelativistic particle of 
mass m = I moving in one dimension in a potential V{x) is given by 

H = y^ + V{x), V{x) = {x'^-3)e-''^^. (2.1) 

This potential has an attractive core surrounded by a repulsive barrier, as shown in Fig. [T[ The 
infinite- volume spectrum of this Hamiltonian for energies £ < 4 is shown in Fig. [T[ The ground 
state is a bound state of even parity, and there is one bound state in the odd parity channel. A 
continuum of scattering states is found for £" > 0, with a narrow resonance in the even-parity 
channel and a broad resonance in the odd-parity sector, both below 4. Even (+) and odd (— ) parity 
scattering phase shifts 5±{E) can be defined in the usual way as the phase between the transmitted 
and incident wave, appearing in the asymptotic wave functions as 

<p^^-'(x) =c+cos|/:|;c|+5+(^) j, (pj^^\x) = c_sgn{x) sin ( k\x\ + 5^ (k)] , (|;c|^oo), (2.2) 
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Figure 1: (Left) The solid curve shows the potential V{x) of the example given in Eq. (2.1 ). The spectrum 
for energies £ < 4 is also shown. The ground state is an even-parity bound state, and there is one bound 
state in the odd-parity channel. Two resonances (a narrow one in the even-parity and a very broad one in the 
odd parity sector) in a continuum of scattering states E > are shown. (Right) The scattering phases 8±{E) 
for this example are shown for E <4. The locations of the resonance energies are shown as dashed vertical 
lines. 



where k = \/2E. These phase shifts are also shown in Fig. [l[ The narrow even-parity resonance is 
seen as a sudden dramatic increase in the phase shift (by an amount comparable to tt) and the broad 
odd-parity resonance appears as a not-so-sudden increase in the phase shift. The resonance masses 
and widths can be extracted by fitting a Breit-Wigner plus a polynomial background to d5±/dE in 
the vicinity of the resonances (or by employing the complex rotation method described in Ref. [[T]]). 

Now consider solving the system in a box of length L such that —^L<x< \L, assuming 
periodic boundary conditions. Note that the potential is now Vi{x) = Y^^^^o^ {x + nL). The infinite 
volume gets tiled into L-length strips in which the potential is replicated. We assume that L is large 
compared to the extent of the potential V so that interactions with mirror potentials is negligible. 

In a finite-box with periodic boundary conditions, the momentum is quantized, so the entire 
spectrum is a series of discrete energies, even for E > 0. The periodic -box spectrum can be de- 
termined in two ways: diagonalization of the Hamiltonian in a basis of states having appropriate 
boundary conditions, and by solving the differential equation and matching to an asymptotic form 
having the correct boundary conditions. Either way, one finds the spectrum shown in Fig. ^ The 
light dotted lines indicate the spectrum for V = which have values In^rp- /L^ for n = 1,2,3, ... , 
plus an « = line in the even parity channel. One sees that a resonance shows up as a series of 
avoided level crossings when viewed against box length L. A narrow resonance, as in the even 
parity sector, can be easily identified by a closely-avoided level crossing, but a broad resonance, as 
in the odd parity sector, is essentially impossible to recognize. 

These plots illustrate the difficulty in extracting resonance parameters from finite-box energies. 
Under certain special circumstances, resonance parameters can be ferreted out using tricks such as 
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Figure 2: (Left) Spectrum of even-parity states for the example Hamiltonian given in Eq. (2.1) in a box of 
extent L with periodic boundary conditions. (Right) Spectrum of odd-parity states. The light dotted lines 
indicate the energy levels for V = which have values 2;r^n^/L^ for n — 1,2,3, . . . , plus an n = line in 
the even parity channel. The horizontal dashed lines show the locations of the infinite-volume resonance 
energies. 



described in Ref. [^. Examining the spectrum in several volumes is important, and knowing the 
pattern of multi-hadron states based on mass determinations of the stable particles and group- 
theoretical combinations of the constituents having total zero-momentum certainly helps. For high 
precision, the phase-shift method of Refs. ^ can be used. In this method, the finite-volume 
energies are used to determine the scattering phase shifts of the partial waves, from which one 
can deduce resonance masses and widths. This method has recently been applied to study the p- 
meson resonance[^]. The kk phase shift was extracted from the finite- volume spectrum, and the 




Figure 3: (Left) mt scattering phase shift showing the p-resonance [|^]. A width Tp = 200^\qq MeV is 
obtained from a Breit-Wigner fit. (Right) Illustration of the binning method of Ref. ^ in which a proba- 
bility distribution is used to identify resonant structure. For various values of L, energies are collected into 
momentum bins to construct a probability distribution W{p). 
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Figure 4: (Left) Results of applying the binning method of Ref. [|6[| near the even-parity resonance of the 
example Hamiltonian in Eq. (2A). Location of the resonance is indicated by the vertical dashed hne. (Right) 
Results for the very broad odd-parity resonance. 



p resonance parameters extracted (see Fig. ^). However, tiie practicality of the pliase-sliift metliod 
has yet to be demonstrated for other resonances in QCD. 

A new histogram method has been recently proposed!^] and applied to synthetic pion-nucleon 
data which mocks up the A-resonance. Energies for several volumes are collected into momen- 
tum bins, and with suitable normalization and free-result subtraction, such histograms produce a 
probability distribution W{p) which shows peaks corresponding to the resonances of interest (see 
Fig. ^). The method has no prior theoretical bias and provides the possibility of seeing resonant 
structure even when avoided level crossings are washed out by a broad resonance. The application 
of this method to the example Hamiltonian in Eq. (2J_) is shown in Fig. Q The narrow resonance 
appears quite clearly, and amazingly, the odd-parity broad resonance is also correctly reproduced. 

Deducing resonance parameters from finite-box spectra remains a difficult challenge, espe- 
cially considering that higher-lying resonances will he above three-particle and four-particle thresh- 
olds and that these resonances can have multiple decay channels. Certainly, further work in this 
area is needed. Perhaps matching the finite-box spectra to that of an effective theory, such as a 
one-boson exchange model, might ultimately be the way to make progress. 



3. Excited stationary states: recent results 

Before discussing the issues and challenges in extracting the energies of stationary states in a 
box, I would hke to summarize the excited-state results which have appeared since the last lattice 
conference. 

A first glimpse of the higher-lying nucleon spectrum in lattice QCD was provided by the 
Hadron Spectrum Collaboration in Ref. [^]. These first results, shown in Fig. were on a small 
12^ X 48 anisotropic quenched lattice with a very heavy pion. Results for both the nucleons and 
A-resonances on 239 quenched configurations on a 16^ x 64 lattice and 167 quenched configura- 
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Figure 5: (Left) Nucleon spectrum from 200 quenched configurations on a 12^ x 48 anisotropic lattice using 
the Wilson gauge and quark actions with Oj ^ 0. 1 fm, Us/at ~ 3.0 and m^i ^ 700 MeV from Ref. ||^]. (Right) 
Nucleon spectrum from 430 Nf ~2 configurations on a 24^ x 64 lattice using a stout-smeared clover fermion 
action and Symanzik-improved gauge action with Us ~ 0.1 fm, a^/at ~ 3, and m^^ — 400 MeV from Ref. [0]. 



tions on a 24^ x 64 lattice using an anisotropic Wilson action with spatial spacing ~ 0.1 fm, 
fli/af ~ 3, and a pion mass niji ~ 490 MeV appeared during the past year[^]. These masses have 
been determined in the past year using 430 Nf = 2 configurations on a 24^ x 64 lattice with a 
stout-smeared clover fermion action and a Symanzik-improved anisotropic gauge action|^. The 
results for a pion mass mj^ = 400 MeV, spacing a., ~ 0. 1 fm and as/ at ~ 3 are shown in Fig. |5[ The 
low-lying odd-parity band shows the exact number of states in each channel as expected from ex- 
periment. The two figures show the splittings in the band increasing as the quark mass is decreased. 
At these heavy pion masses, the first excited state in the G\g channel is significantly higher than 
the experimentally measured Roper resonance. It remains to be seen whether or not this level will 
drop down as the pion mass is further decreased. Most of the levels in the right-hand plot lie very 
close to two-particle thresholds. The use of two-hadron operators will be needed to go to lighter 
pion masses. 

During the past year, extractions of excited meson states have been presented in Ref. 
Results in the pseudoscalar, vector, and axial-vector channels are shown in Fig. |6[ These results 
were obtained using 99 quenched configurations on a 16^ x 32 isotropic lattice with a chirally- 
improved fermion action and the Luscher-Weisz gauge action for lattice spacing a., ~ 0.15 fm and 
a range of pion masses. This work emphasizes the use of derivative sources in correlation matrices 
to obtain the excited states. A search for Ught scalar tetraquark states with isospin 7 = 0,^ was also 
presented at this conference[l 1 ]. 
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Figure 6: (Upper left) Results from Ref. JlO| ] for the first-excited pseudoscalar ^ meson mass against 
pion mass squared from 99 quenched configurations on a 16^ x 32 lattice with a chirally-improved fermion 
action and the Luscher-Weisz gauge action for lattice spacing as 0.15 fm. (Upper right) First-excited 1 
meson mass. (Lower left) Second-excited 1 meson mass. (Lower right) Ground and first-excited 1++ 
meson masses. 

4. Excited stationary states: key issues 

Reliably capturing the masses of excited states requires the computation of correlation matrices 
Cij{t) = (0|r<I>;(?)<I>j(0)|0) associated with a large set of N different operators <I>,(0- It has been 
shown in Ref. [jl^ that the N principal effective masses Wa{t), defined by 

where Xait,to) are the eigenvalues of C(fo)^''^ C{t) C(fo)^^^ and to < t/2 is usually chosen, 
tend to the eigenenergies of the lowest A'^ states with which the N operators overlap as t becomes 
large. The eigenvectors associated with Xa{t,to) can be viewed as variationally optimized opera- 
tors. When combined with appropriate fitting and analysis methods, such variational techniques 
are a particularly powerful tool for investigating excitation spectra. To extract the stationary state 
energies, one can fit a single-exponential or a sum of two exponentials to each principal correla- 
tor; alternatively, optimized operators can be determined on an early time slice, and a fit to the 
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correlation matrix of optimized operators can be carried out. The use of both methods is a good 
consistency check. 

The use of operators whose correlation functions C{t) attain their asymptotic form as quickly 
as possible is crucial for reliably extracting excited hadron masses. An important ingredient in con- 
structing such hadron operators is the use of smeared fields. Operators constructed from smeared 
fields have dramatically reduced mixings with the high frequency modes of the theory. Both link- 
smearing and quark-field smearing should be applied. Since excited hadrons are expected to be 
large objects, the use of spatially extended operators is another key ingredient in the operator de- 



sign and implementation. A more detailed discussion of these issues can be found in Ref. [13] 



Spatial links can be smeared using the stout-link procedure described in Ref. [14]. The stout- 
link smearing scheme is analytic, efficient, and produces smeared links which are automatically 
elements of SU (3) without the need for a projection back into SU (3). Note that only spatial staples 
are used in the Unk smoothening; no temporal staples are used, and the temporal link variables are 
not smeared. The smeared quark fields can be defined by 



V{x)=^l + ^Aj (4.1) 

where and ria are tunable parameters {ria is a positive integer) and the three-dimensional covari- 
ant Laplacian operators are defined in terms of the smeared link variables Uj{x) as follows: 

AO{x) = £ (Uk{x)0{x+k) - 0{x)Y (4.2) 

/t=±l,±2,±3V / 

where 0{x) is an operator defined at lattice site x with appropriate color structure, and noting that 
U-k{x) = U^{x—k). The smeared fields \j/ and \j7 are Grassmann- valued; in particular, these fields 
anticommute in the same way that the original fields do, and the square of each smeared field 
vanishes. 

Hadron states are identified by their momentum p, intrinsic spin J, projection A of this spin 
onto some axis, parity P = ±1, and quark flavor content (isospin, strangeness, etc.). Some mesons 
also include G-parity as an identifying quantum number. If one is interested only in the masses of 
these states, one can restrict attention to the p = sector, so operators must be invariant under all 
spatial translations allowed on a cubic lattice. The little group of all symmetry transformations on 
a cubic lattice which leave p = invariant is the octahedral point group Oh, so operators may be 
classified using the irreducible representations (irreps) of Oh. For mesons, there are ten irreducible 
representations A\g,A2g,Eg,T\g, T2g,Aui,A2u,Eu,T\u,T2u- The representations with a subscript g{u) 
are even (odd) under parity. The A irreps are one dimensional, the E irreps are two dimensional, and 
the T irreps are three-dimensional. The Ai irreps contain the 7 = 0,4,6, 8, . . . states, the A2 irreps 
contain the 7 = 3, 6, 7, 9,... states, the E irreps contain the 7 = 2, 4, 5, 6, 7,... states, the T\ irreps 
contain the spin 7 = 1,3,4,5, . . . mesons, and the T2 irreps contain the spin 7 = 2,3,4,5, . . . states. 
For baryons, there are four two-dimensional irreps G\g,G\u,G2g, G2u and two four-dimensional 
representations Hg and //„. The G\ irrep contains the J = j,^,^,^ , ■ ■ ■ states, the H irrep contains 
the 7 = - states, and the G2 irrep contains the 7 = .. states. The continuum- 

Umit spins 7 of our states must be deduced by examining degeneracy patterns across the different 
Oh irreps. 
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Figure 7: The spatial arrangements of the extended three-quark baryon operators. Smeared quark-fields 
are shown by solid circles, line segments indicate gauge-covariant displacements, and each hollow circle 
indicates the location of a Levi-Civita color coupling. For simplicity, all displacements have the same length 
in an operator Results presented here used displacement lengths of Saj 0.3 fm). 



The authors of Ref . JTsI ] advocate operators designed with one eye towards maximizing over- 
laps with the low-lying states of interest, and the other eye towards minimizing the number of 
sources needed to calculate the required quark propagators. They emphasize that a construction 
method which can be easily adapted for baryons, mesons, hybrid states, and multi-hadron systems 
is ideal. Since the calculation of quark propagators can be computationally expensive, baryon, 
meson, and multiquark operators which share the same basic building blocks is recommended. 

Thus, these authors advocate a three-stage approach to constructing hadron operators. First, 
basic building blocks are chosen. These are taken to be smeared covariantly-displaced quark fields 

iD?v)lAWDfraa^ -3<j<3, (4.3) 

where A is a flavor index, a is a color index, a is a Dirac spin index, and the /^-link gauge-covariant 
displacement operator in the j-th direction is defined by 

D'^f\x,x') = Uj{x) Uj{x+j)...Uj{x+{p-l)j)5,,^^^^,j, b[P\x,x') = 5,,,, (4.4) 

for 7 = ±1 , ±2, ±3 and p>\, and where j = defines a zero-displacement operator to indicate no 
displacement. Next, elemental operators Bf{t,x) are devised having the appropriate flavor structure 
characterized by isospin, strangeness, etc., and color structure constrained by gauge invariance. For 
zero momentum states, translational invariance is imposed: Bf{t) = T^x^f (fi^)- Finally, group- 
theoretical projections are applied to obtain operators which transform irreducibly under all lattice 
rotation and reflection symmetries: 

^AAF(,) = ^ I r^tliR) Ur Bfit) Ul (4.5) 



So? 



r, ReOf, 



where is the double group of Oh, R denotes an element of Of^, gpo is the number of elements in 
O^, and dx is the dimension of the A irreducible representation. Projections onto both the single- 
valued and double- valued irreps of Oh require using the double group in Eq. (4.5). Given Mb 



0. - U U R 

single- singly- doubly- triply- triply- 

site displaced displaced-L displaced-U displaced-0 

Figure 8: The spatial arrangements of the quark-antiquark meson operators. In the illustrations, the smeared 

quarks fields are depicted by solid circles, each hollow circle indicates a smeared "barred" antiquark field, 

and the solid line segments indicate covariant displacements. 
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Figure 9: Effective masses M{t) for unsmeared (black circles) and smeared (red triangles) operators 
Oss, OsD, Otdt, which are representative single-site, singly-displaced, and triply-displaced-T nucleon 
operators, respectively. Top row: only quark-field smearing = 32, Os ~ 4.0 is used. Middle row: 
only link-variable smearing np — 16, rtpp — 2.5 is applied. Bottom row: both quark and link smearing 
«(j = 32, (7v = 4.0, Hp = 16, «pP = 2.5 are used, dramatically improving the signal for all three operators. 
Results are based on 50 quenched configurations on a 12^ x 48 anisotropic lattice using the Wilson action 
with Us ~ 0.1 fm, Qs/at ~ 3.0. 



elemental Bf operators, many of the projections in Eq. ( p~5| ) vanish or lead to linearly-dependent 
operators, so one must then choose suitable linear combinations of the projected operators to obtain 
a final set of independent baryon operators. Thus, in each symmetry channel, one ends up with a 
set of r operators given in terms of a linear superposition of the Mb elemental operators. The 
different spatial configurations (see Fig. ^ for the baryon configurations and Fig. |8| for the meson 
configurations) yield operators which effectively build up the necessary orbital and radial structures 
of the hadron excitations. The design of these operators is such that a large number of them can 
be evaluated very efficiently, and components in their construction can be used for both meson, 
baryon, and multi-hadron computations. 

Finding appropriate smearing parameters is a first crucial part of any hadron spectrum calcu- 
lation. Fig. ^ demonstrates that both quark-field and link-field smearing are needed in order for 
spatially-extended baryon operators to be useful[15]. It is important to use the smeared links when 
smearing the quark field. Link smearing dramatically reduces the statistical errors in the correla- 
tors of the displaced operators, while quark-field smearing dramatically reduces the excited-state 
contamination. 

The above approach to designing hadron and multi-hadron interpolating fields leads to a very 
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Figure 10: Effective masses for the lowest eight levels in the Gig channel (left) and //„ channel (right). 
The green circles are fixed-coefficient effective masses, whereas the red squares are the principal effective 
masses. Fit values are shown by the blue lines. Results are based on 200 quenched configurations on a 
12^ X 48 anisotropic lattice using the Wilson gauge and quark actions with Oj ^0.1 fm, Us/at ^ 3.0 and 
nij, - 700 MeV. 



large number of operators. It is not feasible to do spectrum computations using all of the operators 
so designed; for example, in the Gig symmetry channel for nucleons, the above procedure leads to 
179 operators. It is necessary to prune down the number of operators. Six months of exploratory 
testing and trials led to the following guideline: noise is the enemy, so a procedure that keeps a 
variety of operators while minimizing the effects of noise works best. Some operators are intrin- 
sically noisy and must be removed. In addition, a set of operators, each with little intrinsic noise, 
can allow noise to creep in if they are not sufficiently independent of one another. 

In Ref. [^, the following procedure is advocated. (1) First, remove operators with excessive 
intrinsic noise. This can be done by examining the diagonal elements of the correlation matrix and 
discarding those operators whose self-correlators have relative errors above some threshold for a 
range of temporal separations. Of course, this requires a low-statistics Monte Carlo computation on 
a reasonably small lattice. (2) Second, prune within operator types (single-site, singly-displaced. 
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etc.) based on the condition number of the submatrices 

Ciiit) = , ^'^^^^ :, t = at. 

The condition number is taken to be the ratio of the largest eigenvalue over the smallest eigenvalue. 
A value near unity is ideal. For each operator type, the set of about six operators which yields the 
lowest condition number of the above submatrix is retained. (3) Lastly, prune across all operator 
types based again on the condition number of the remaining submatrix as defined above. In this 
last step, the goal is to choose about 16 operators, keeping two or three of each type, such that a 
condition number reasonably close to unity is obtained. As long as a good variety of operators is 
retained, the resulting spectrum seems to be fairly independent of the exact choice of operators at 
this stage. Eigenvectors from a variational study of the operators can also be used to fine tune the 
choice of operators. 

Calculations in Ref. using about 16 operators in all irreps for the nucleon channel on 
only 200 configurations were very successful. The extraction of 8 energy levels in each irrep was 
possible, which was a major milestone achieved. A determination of the hadron spectrum requires 
the ability to extract several excited energy levels, and up until that time, it was not known whether 



or not extracting more than one or two levels would be possible. Fig. |10| shows the signal quality 
in the Gig and //„ irreps for the nucleon excitations from that first calculation. Similar calculations 
for the A resonance spectra have also been achieved[|T6|, [l7| ]. Comparison of these results with 
experiment is not justified since the quenched approximation was used, an unphysically large u,d 
quark mass was used, and the lattice volume is too small. 

It is my strong opinion that the use of correlation matrices is the best way to extract excited- 
state energies reliably. However, there are efforts to deduce information about excited states from 
single correlation functions. Bayesian statistics have been used[18, [T^, as well as maximum en- 



tropy methods 1 2C, 21]. A novel evolutionary fitting method has been proposed[E3], and a new 



method based on statistical concepts which relies heavily on simulation techniques was presented 



at this conference! 23 ] . 



5. Stochastic estimates of many-to-many quark propagators with source dilution 
variance reduction 

To study a particular eigenstate of interest, all eigenstates lying below that state must first be 
extracted, and as the pion gets lighter in lattice QCD simulations, more and more multi-hadron 
states will lie below the excited resonances. Consider a baryon at rest. An appropriate quantum 
operator for a baryon at rest typically has the form 

B{p = 0,t) = ^'£(pB{r,t), (5.1) 

' r 

where V is the volume of the lattice and <pfi(r,?) is an appropriate localized interpolating field. In 
the above equation, the summation over spatial lattices makes the operator translationally invariant, 
producing a zero momentum state. A baryon correlator, thus, has a double summation over spatial 
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sites: 

(0|B(p = 0,0B(P = 0,0)10) = -^£(0|<Pb(^,0^b(>',0)|0). (5.2) 

x,y 

Evaluating the above correlator requires computing the y —^x element of the quark propagators. In 
other words, the quark propagators from all spatial sites y on time slice f = to all spatial sites x on 
later time slice t >0 must be known. Computing all such elements of the propagators exactly is not 
possible (except on very small lattices). In the example above, this problem can be circumvented 
by appealing to translational invariance to limit the summation over the source site to a single site: 

{0\B{p = 0,t)B{p = 0,0)|0) = ^£(0|(pB(a:,09B(0,0)|0). (5.3) 
However, a good baryon-meson operator of total zero momentum typically has the form 

B{p,t)M{-p,t) = ^Y.'PBM'PMiy^ty'"^''''^ (5-4) 

where (pM{y,t) is a localized interpolating field for a meson. In the evaluation of the temporal cor- 
relations of such a multi-hadron operator, it is not possible to completely remove all summations 
over the source site. Hence, the need for estimates of the quark propagators from all spatial sites on 
a time slice to all spatial sites on another time slice cannot be sidestepped. Ultimately, some cor- 
relators will involve disconnected diagrams which necessarily involve all-to-all quark propagators. 
Hence, all-to-all (or many-to-many) quark propagators are becoming mandatory, and some way of 
stochastically estimating them is needed. 

Random noise vectors Tj whose expectations satisfy E{rii) = and ^(TjiTjp = 5ij are useful 
for stochastically estimating the inverse of a large matrix M as follows. Assume that for each of Nr 
noise vectors, we can solve the following linear system of equations: MX''' = T]^'') for X^''\ Then 
ZW =M-iT]W,and 

£(X,-T];) =£(£m,^1t],t];) =£M,^i£(T],T];) =Xm,^^4; =M,^i. (5.5) 

k k k 

The expectation value on the left-hand can be approximated using the Monte Carlo method. Hence, 
a Monte Carlo estimate of M^^ is given by 

Mr^ ^ Urn —Vxy'r^)'''*, where MXW = t]W. (5.6) 

Unfortunately, this equation usually produces stochastic estimates with variances which are much 
too large to be useful. 

Progress is only possible if stochastic estimates of the quark propagators with reduced vari- 
ances can be made. Techniques of diluting the noise vectors have been developed which accomplish 
such a variance reduction! 24 , 25, 26, 27, 28, 29]. A given dilution scheme can be viewed as the 
application of a complete set of projection operators. To see how dilution works, consider a gen- 
eral N X N matrix M having matrix elements Mij. Define some complete set of N xN projection 
matrices P^"^ which satisfy 

p[a) p(h) _ ^ab p{a) ^ pia)"! _ p{a) ^ ^^j-^ 
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Then observe that 

a a a 

= WP^Eir^^^})Pfl = We (pSn.n}pfl) • (5.8) 

a a 

Define 

[a]_p(a) *p(fl)_p(a)* * 

and further define X^"^ as the solution of 

M,,xi"l=77fl, (5.10) 

then we have 

a a 

Although the expected value of Y^a'Hk'^'nf^* the same as TjjtTj*, the variance of Y.a'Hk'^'nf^* 
significantly smaller than that of rjkrij. For both Z4 and U{\) noise, we have 

Vflr(Re(77,-T];)) = Vflr(lm(T],-77;)) = ^(1 - 5,-y). 

Although the variance is zero for / = j, there is a significant variance for all / 7^ j. The dilution 
projections ensure exact zeros for many of the off-diagonal elements, instead of values that are only 
statistically zero. In other words, many of the / 7^ j elements become exactly zero. 

Of course, the effectiveness of the variance reduction depends on the projectors chosen. A 
particularly important dilution scheme for measuring temporal correlations in hadronic quantities 
is "time dilution" where the noise vector is broken up into pieces which only have support on a 
single time slice: 

Paa;bp(.^^^'^y^^') = ^"bSai}Sxy5Bt8Bt', B = 0, 1 , • • • ,A^, - 1 , (time dilution) , (5.12) 

where Nf is the number of time slices on the lattice, a,b are color indices, and a, jS are spin indices. 
Spin and color dilution are two other easy-to-implement schemes: 

Paa-bpi^^^'y^^') = ^"b^BaSsiiSxydtt', 5 = 0,1,2,3, (spin dilution), (5.13) 
Paa;bpi^^f'^y^f') = ^BaSBbSapSxyStt', B = 0,1,2, (color dilution). (5.14) 

Various spatial dilution schemes are possible, too. For example, even-odd dilutions are simple to 
implement. The above dilution projectors can also be combined to make hybrid schemes. 

Before presenting tests of these different dilution schemes, an important remark about the use 
of stochastic quark propagators should be mentioned. The use of Eq. ( |5.11| ) to approximate quark 
propagators leads to a very desirable source-sink factorization. Consider a baryon correlator of the 
form 

C„ = 4<$Qli>Q^d£. (5.15) 
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Figure 11: (Upper row) Effective masses for a single-site (left), singly-displaced (middle), and triply- 
displaced-T (right) nucleon operator using quark propagators evaluated with the standard point-to-all 
method. (Lower row) The effective masses for the same nucleon operators but using stochastic quark propa- 
gators with time-Hspin-Hcolor dilution. Without dilutions, the errors in these effective masses would be orders 
of magnitude larger. An effective mass defined using a time separation 3fl, is used in these plots. These re- 
sults used 100 quenched configurations on an anisotropic 12^ x 48 lattice with a Wilson fermion and gauge 
action. 



where Q^^^ denotes a quark propagator of flavor A and all other quark indices have been combined 
into a single index / or j, and so on. Stochastic estimates of this correlator using Eq. ( 5.11 ) lead to 
the form 

1 



d All Bile 



ri- 



{Br)[ds]^{Br)[dB]* 
j 



9j -'m 



JCr)[dc] {Cr)[dc]* 



(5.16) 



where r labels the noise vectors, dA,dB,dc are the dilution indices, rj are the noise vectors, and (p 
are the solution vectors. If one defines 



.{r)ldAdBdc] _ U).JAr)[dA] {Br)[dB] (Cr)[dc] 



{r)[dAdBdc] 



(5.17) 

"/ — ^iik'ii 'li 'Ik ' (5.18) 

then the baryon correlator becomes a glorified dot product of the source vector with the sink vector: 



M „{Arm (Br)ldB] {Cr)ldc] 
^ijk'li 'Ij 'Ik ' 



1 



r d^dBdc 



(r) \dAdBdc\^r) [dAdgdc] * 



(5.19) 



The source and sink vectors in Eqs. (5.17) and (5.18) can be separately evaluated for a variety of 
operators, and the dot product applied afterwards to evaluate the matrix of correlation functions. 
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Figure 12: (Left) The relative errors in the correlation function of a single-site nucleon operator for temporal 
separation t = 5at evaluated using stochastically-estimated quark propagators with different dilution schemes 

1 /2 

against 1 /N^^^ , where A^inv is the number of Dirac matrix inversions required. The open circle shows the 
point-to-all error, and the horizontal dashed line shows the gauge-noise limit. The black (red) dashed-dotted 
line shows the decrease in error expected by simply increasing the number of noise vectors, starting from 
the time (time + even/odd-space) dilution point. (Right) Same as the left plot, except for a triply-displaced-T 
nucleon operator. These results used 100 quenched configurations on an anisotropic 12^ x 48 lattice with a 
Wilson fermion and gauge action. 



Different ABC permutations of the noise vectors must be stored in order to accommodate all needed 
Wick contractions. The use of stochastic all-to-all quark propagators has led to an enormous simpli- 
fication of the effort required to compute the hadron correlation matrices through this source-sink 
factorization. Another advantage of this approach is the fact that, given suitable non-zero momenta, 
these same baryon and meson operators can be combined later to make multi-hadron operators. 

The effectiveness of stochastically-estimated all-to-all quark propagators using diluted noise 
vectors is demonstrated in Fig. [ri[ This figure compares the effective masses for a single-site, 
singly-displaced, and triply-displaced-T nucleon operator using quark propagators evaluated with 
the conventional point-to-all method (top row) and with the all-to-all stochastic method including 
time-i-spin-i-color dilutions (bottom row). The fact that these effective masses have comparable er- 
rors indicates that the stochastic method with suitable dilutions has not introduced any appreciable 
noise into the final mass extractions. 



A comparison of different dilution schemes has been presented at this conference |gOp. Fig. |12 
shows the relative errors in the correlation functions of a single-site and a triply-displaced-T nu- 
cleon operator for temporal separation t = 5at evaluated using stochastically-estimated quark prop- 
agators with different dilution schemes against 1 /N^^^ , where A^i^y is the number of matrix inver- 
sions required. These results were obtained using 100 quenched configurations on an anisotropic 
12^ X 48 lattice with a Wilson fermion and gauge action. The open circles show the point-to-all 
errors, and the horizontal dashed lines show the gauge-noise limits. The black (red) dashed-dotted 
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lines show the decrease in error expected by simply increasing the number of noise vectors, start- 
ing from the time (time + even/odd-space) dilution points. These computations are dominated by 
the inversions of the Dirac matrix, so using the number of matrix inversions A^inv to compare com- 
putational efforts is reasonably fair. The advantage in using increased dilutions compared to an 
increased number of noise vectors with only time dilution is evident in the plots. However, this 
advantage quickly diminishes after time -i- even/odd-space dilution, or time-i-color, or time-i-spin 
dilution. Note that time-i-spin-i-color-i-even/odd-space dilution yields an error comparable with the 
gauge-noise limit using only a single noise vector! First results for multi-hadron operators were 



also presented at this conference! 31 



These encouraging results demonstrate that the inclusion of good multi-hadron operators will 
certainly be possible using stochastic all-to-all quark propagators with diluted-source variance re- 
duction. In fact, just before this conference, the authors of Ref. began exploring a new method 
that might allow nearly-exact determinations of many-to-many quark propagators without intro- 
ducing any noise vectors at all. The method exploits a novel, cleverly-devised choice of quark-field 
smearing to facilitate the nearly-exact computations. Details and tests of this method should appear 
very soon. 



6. Summary and outlook 

This talk discussed the key issues and challenges in exploring excited hadrons in lattice QCD. 
The importance of multi-hadron operators and the need for all-to-all quark propagators were em- 
phasized. The challenge of dealing with unstable states (resonances) in a box was outlined, and the 
technology associated with extracting excited stationary-state energies, including operator design 
and field smearing, was detailed. Efforts in variance reduction of stochastically-estimated all-to- 
all quark propagators using source dilutions were described, and results on excited hadrons which 
appeared during the last year were summarized. 

Given the major experimental efforts to map out the QCD resonance spectrum, such as Hall 
B and the proposed Hall D at Jefferson Lab, ELSA associated with the University of Bonn, COM- 
PASS at CERN, PANDA at GSl, and BESlll in Beijing, there is a great need for ab initio deter- 
minations of such states in lattice QCD. The exploration of excited hadrons in lattice QCD is well 
underway. 

This work was supported by the National Science Foundation through awards PHY 0653315 
and PHY 0510020. 
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